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A  RAPID  METHOD  FOR  ESTIMATING  EQUIVALENT  PARABOLIC  POLARS  FROM 

WIND-TUNNEL  TESTS 

by 

D.  E.  Lean 
S.  P.  Fiddes 


SUMMARY 


A  method  for  obtaining  a  'best-fit*  parabolic  drag  polar  from  measured 
force  results  is  presented.  The  method  is  used  in  estimating  the  blockage 
correction  due  to  separated  flow  that  must  be  applied  to  wind-tunnel  results. 

An  important  part  of  the  method  is  the  use  of  computer  graphics  to  identify 
the  part  of  the  measured  drag  polar  that  may  best  be  described  by  a  parabola, 
and  thus  used  as  the  basis  for  a  curve  fit. 


The  method  is  fast,  when  compared  to  existing  curve-fitting  techniques, 
and  suited  to  use  with  an  'online'  data  reduction  system. 
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INTRODUCTION 


In  the  wind-tunnel  testing  of  lifting  models,  it  is  often  necessary  to  estimate  the 
profile  and  induced  drag  components  of  the  total  drag,  the  remainder  being  referred  to  as 
the  drag  due  to  separated  flow.  A  particular  need  for  this  form  of  drag  decomposition 
arises  in  the  application  of  constraint  corrections  to  wind-tunnel  results,  in  particular 
the  correction  for  separated  wake  blockage  (see,  for  example.  Ref  I).  The  presence  of 
a  large  amount  of  separated  flow  over  the  model,  inferred  from  non-zero  values  of  drag 
due  to  separated  flow,  leads  to  the  generation  of  a  separated  flow  wake.  The  constraining 
effect  of  the  tunnel  walls  means  that  the  separated  flow  wake  can  only  be  accommodated  in 
the  tunnel  by  an  acceleration  of  the  flow  in  the  vicinity  of  the  model.  This  acceleration 
of  the  flow  is  regarded  as  a  change  in  freestream  reference  kinetic  pressure  (Aq)  which 
must  be  found  if  the  aerodynamic  coefficients  are  to  be  calculated  correctly. 

To  make  the  drag  decomposition  a  more  concrete  process,  it  is  assumed  that  in 
attached  flow  the  drag  of  a  general  wing  (or  wing-body  combination)  may  be  written  as 


(1) 


where  is  the  'attached  flow  drag  coefficient', 

A 

k  is  the  induced  drag  factor, 

AR  is  the  aspect  ratio  of  the  model  wing 

and  C.  ,  C_  are  constants. 

The  parabolic  drag  polar  given  by  equation  (t)  is  an  idealisation,  having  a  rigorous 

foundation  only  in  linear  finite-wing  theory,  but  it  has  been  found  to  be  a  useful  approxi¬ 

mation  for  the  drag  of  a  model  without  extensive  separation.  However,  in  these  cases, 
the  lift-dependent  drag  term  includes  the  variation  of  profile  drag  with  sectional  lift 
coefficient.  This  profile  drag  variation  should,  in  principle,  be  extracted  so  that  the 
blockage  correction  due  to  the  'attached  flow  wake'  may  be  estimated,  leaving  the  remain¬ 
ing  lift-dependent  drag  as  the  induced  drag  alone. 

At  large  lifts  (high  C^s)  the  measured  drag  (C^)  will  generally  exceed  the  'attached 
flow'  drag  given  by  ,  and  this  difference  is  attributed  to  the  separated  flow  drag, 

which  is  assigned  a  coefficient  value  of  “  ^j)  “  ^D^  ^^8  *)• 
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Maskell  has  given  an  expression  relating  the  separated  flow  drag  coefficient  and  the 
blockage  it  induces,  via: 


q 


where  Aq/q  is  the  wake  blockage  factor  used  to  correct  the  tunnel  reference  kinetic  press¬ 
ure,  S  is  the  reference  area  of  the  model,  and  C  is  the  cross-sectional  area  of  the 
tunnel  working  section. 


The  problem  addressed  in  this  paper  is  the  estimation  of  ,  or,  more  precisely, 

the  estimation  of  the  constants  Cqqi  and  k  in  equation  (1),  from  experimental 

data,  so  that  the  separated  flow  drag  coefficient  may  be  found.  It  must  be  emphasised,  in 
view  of  the  remarks  above  concerning  the  variation  of  profile  drag  coefficient  with  lift, 
that  the  resulting  values  for  k  do  not  represent  the  induced  drag  factor.  The  results 
from  this  analysis  should  not,  therefore,  be  regarded  as  indicating  the  lifting  efficiency 
of  a  particular  wing. 

2  PREVIOUS  METHODS 

Inspection  of  equation  (1)  reveals  that  the  minimum  value  of  Cq^  is  Cq^  ,  and  it 
thus  appears  that  to  determine  Cq^  all  that  is  required  is  to  determine  the  minimum  drag 
measured  in  the  experiment.  However,  this  assumes  that  the  observed  minimum  drag  lies  in 
the  range  of  lift  for  which  equation  (1)  is  an  adequate  description  of  the  observed  drag 
polar.  This  is  not  generally  the  case,  as  for  some  highly-cambered  wings  -  such  as  a  high- 
lift  wing  -  minimum  drag  can  occur  at  relatively  high  lift  coefficients,  as  separations 
appear  on  the  lower  surface  of  a  slat,  leading  to  a  drag  rise  as  the  lift  decreases.  An 
example  of  such  behaviour  will  be  given  later.  This  highlights  a  particular  difficulty 
with  the  problem  of  estimating  an  'equivalent'  parabolic  drag  polar  -  it  does  not  apply  to 
all  of  the  observed  drag  polar,  and  the  region  where  it  is  applicable  must  be  estimated  in 
some  way.  Early  attempts  at  fitting  a  parabolic  polar  to  part  of  the  measured  polar  relied 
on  a  trial  and  error  procedure  of  the  following  type.  First,  an  estimate  is  made  of  Ci^  , 
and  a  curve  of  against  -  CLq)  is  drawn.  If  the  guessed  value  of  is  close 

to  the  actual  value,  then  the  curve  should  be  a  substantially  straight  line.  A  rather 
arbitrary  decision  is  made  on  what  constitutes  a  satisfactory  'degree  of  straightness', 
and  the  process  repeated  with  a  new  guess  for  until  a  satisfactory  fit  is  obtained. 

From  the  slope  and  intercept  of  this  curve,  values  for  k  and  Cqq  may  be  obtained. 
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Smith  attempts  to  rationalise  the  process  by  means  of  a  sequence  of  least-squares 
curve  fits.  Again,  a  range  of  values  is  estimated,  and  for  each  a  least- 

squares  fit  is  taken  over  a  defined  range  of  lift  coefficients  on  the  observed  polar. 

A  correlation  coefficient  is  then  calculated  which  gives  a  measure  of  the  goodness  of  fit 
of  the  curve  to  the  specified  part  of  the  polar.  is  then  varied  over  the  range  speci¬ 

fied  in  an  attempt  to  find  a  tna-riimim  value  for  the  correlation  coefficient. 

The  virtue  of  this  method  is  that  once  a  range  of  likely  values  has  been 

supplied,  and  the  part  of  the  measured  drag  polar  thought  to  be  near-parabolic  has  been 
estimated,  the  process  is  well  defined  and  may  be  automated.  However,  some  skill  is 
required  in  estimating  the  parabolic  part  of  the  observed  drag  polar,  for  which  the  method 
gives  no  guidance.  A  similar  problem  arises  in  estimating  a  range  of  values  to  be 

tried.  Furthermore,  in  searching  for  a  maximum  value  of  the  correlation  coefficient, 
a  number  of  curve  fits  has  to  be  taken,  which  may  lead  to  unacceptable  computation  times 
if  a  large  number  of  runs  has  to  be  processed. 

It  was  therefore  felt  desirable  to  develop  a  method  which  relied  less  on  inspired 
guesswork  by  the  user,  and  less  on  trial  and  error  curve  fitting. 
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NEW  METHOD 


Referring  to  equation  Cl)>  consider  two  adjacent  datapoints  (Cl^,  Cq^)  and 

If  these  datapoints  lie  on  a  parabolic  part  of  the  observed  drag  polar, 

we  may  write: 


(2) 


Therefore,  if  we  plot  against  ,  we  will 

obtain  a  straight  line.  Mote  that  because  a  parabolic  drag  polar  is  assumed,  the 
'finite-difference'  expression  is  in  fact  the  derivative 

of  the  drag  polar  at  a  lift  coefficient  of  +■  ,  and  is  thus  written  as 

dC^/dC^  .  In  practice,  we  expect  only  part  of  the  measured  drag  polar  to  be  adequately 
represented  by  a  parabola,  so  only  part  of  the  plot  obtained  in  this  way  will  consist  of 
a  straight  line.  However,  with  the  naked  eye,  it  is  easier  to  identify  a  straight  line 
than  a  parabola,  and  so  the  plot  quickly  iiuiicates  that  part  of  the  measured  drag  polar 
which  is  parabolic,  and  thus  the  part  of  the  polar  that  may  be  used  in  a  curve-fit  to 


obtain  C 


Lq 


etc. 


Furthermore,  the  plot  like  that  described  by  equation  (2)  may  be  used  directly  to 
determine  k  and  ,  as  the  straight-line  portion  of  the  plot  has  a  slope  of  ZR/ttAR  , 

and  an  intercept  at  .  Thus  one  straight-line  curve  fit  (least-squares  if  necessary) 

over  a  range  of  datapoints  that  are  (hopefully)  obvious  from  the  plot  is  sufficient  to 
determine  and  k  .  To  determine  Cq^  ,  the  value  of  Cx,q  found  from  the  plot 

described  above  is  used  to  obtain  a  plot  of  against  (C^  -  CLq)^  •  A  straight-line 

curve  fit  to  this  plot,  over  the  same  range  of  s  used  in  the  first  fit,  will  determine 
Cqq  from  the  intercept  of  the  fitted  curve  with  the  axis. 


The  new  method  is  regarded  as  having  two  distinct  advantages  over  the  previous 
methods: 


(i)  the  range  of  measured  datapoints  over  which  the  parabolic  drag  polar  should  be 
fitted  is  obvious  from  inspection  of  the  first  plot; 

Cii)  no  iterative  curve  fitting  is  required.  Only  two  curve  fits  are  required 
per  case. 


An  internal  consistency  check  is  possible,  as  two  independent  estimates  of  k  are 
available  from  the  slopes  of  the  first  and  second  plots. 
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4  RESULTS 

Results  are  presented  here  using  data  obtained  from  tests  on  a  civil  transport  model 
in  the  RAE  5  metre  Tunnel.  Four  configurations  of  the  model  are  considered,  corresponding 
to  cruise,  take-off  with  two  flap  settings,  and  a  landing  case. 

The  associated  graphical  results  obtained  by  using  the  present  method  are  shown  in 
Figs  2  to  5.  Each  Figure  has  three  parts.  The  upper  part  shows  the  plot  of  dC^/dC^ 
against  that  is  used  to  identify  that  parabolic  part  of  the  polar,  ie  the  part  of 

this  plot  that  is  straight,  or  very  nearly  so.  This  is  easily  seen  for  the  cruise  case 
(Fig  2),  while  the  take-off  and  landing  cases  (Figs  3  to  5)  show  varying  degrees  of 
'straightness'  for  this  curve.  However,  they  all  show  a  characteristic  feature,  and  that 
is  a  straight  region  in  the  middle  of  the  curve  fdiich  has  the  lowest  slope  of  any  part  of 
the  curve.  This  is  the  part  of  the  polar  where  drag  changes  least  rapidly  with  lift,  and 
is  identified  with  the  'attached  flow*  part  of  the  polar.  Either  side  of  this  part  of  the 
polar,  the  drag  rises  rapidly  with  increasing  lift  (on  the  right-hand  side)  or  decreasing 
lift  (on  the  left-hand  side) .  A  straight-line  fit  is  taken  through  the  points  lying  on 
this  middle  portion  of  the  cxirve,  between  the  region  delimited  by  the  vertical  bars.  From 

the  slope  and  intercept  of  the  straight  line,  an  estimate  is  obtained  for  k  and  . 

“  2 

The  value  of  is  then  used  to  produce  the  second  curve,  of  against  (C^  -  CLq) 

This  produces  another  straight  line  for  part  of  the  curve,  and  a  least-squares  straight 
line  fit  is  taken  through  the  same  range  of  lift  coefficients  used  for  the  straight-line 
fit  in  the  first  plot.  The  intercept  of  the  second  straight  line  determines  a  value  for 
Cqq  .  The  parabolic  polar  is  thus  determined.  The  lower  part  of  the  Figure  shows  a  plot 
of  the  original  data  with  the  fitted  polar  also  drawn.  Note  that  Figs  3  to  5  have  C^q 
outside  the  range  of  the  parabolic  part  of  the  curve.  This  is  because  of  the  drag  rise 
with  decreasing  lift  associated  with  separated  flow  on  the  slat,  as  described  earlier. 

3 

These  results  are  compared  to  those  obtained  by  using  Smith's  method  in  Table  I. 

It  is  seen  that  there  is  quite  good  agreement  for  the  values  of  Cqq  ,  but  that  there  are 
marked  differences  in  the  values  of  and  k  for  some  cases.  However,  it  is  mis¬ 
leading  just  to  compare  the  values  of  CjOg  and  k  -  we  are  concerned  with  the 

differences  that  this  implies  in  the  final  wake  blockage  correction.  Table  I  therefore 
shows  the  implied  value  of  the  'separated  flow'  drag  coefficient,  and  the  corresponding 
blockage  correction.  The  final  values  of  wake  blockage  correction  at  obtained  by 

each  method  are  in  good  agreement,  despite  the  differences  in  k  and  used  in  each 

case.  This  is  because  a  'low'  value  of  Cx,^  is  offset  by  a  'low'  value  of  k  .  A  similar 

4  ^ 

result  was  noted  by  Lovell  .  This  indicates  that  the  problem  of  determining  k  and 
is  to  some  degree  ill-conditioned,  in  that  appreciable  differences  in  k  and  Cx.^  do  not 
affect  the  'goodness  of  fit'  very  strongly.  However,  the  converse  of  this  argument  is  that 
the  final  parabolic  polar  fit  is  not  very  sensitive  to  the  actual  values  of  k  and  Cx.^ 
obtained  by  the  curve  fit. 

Therefore;  while  the  estimation  of  accurate  values  of  k  and  Cx.^  from  experimental 
data  is  difficult,  the  estimation  of  the  corresponding  parabolic  polar  is  easier. 
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Two  further  applications  of  the  present  method  are  shown  in  Figs  6  and  7.  Fig  6 
shows  the  analysis  of  some  results  from  testing  a  trainer  model,  while  Fig  7  is  based  on 
results  obtained  from  tests  on  a  canard-delta  combat  aircraft  model.  Both  Figures  show 
an  easily  identifiable  linear  portion  of  the  top  curve  (dCj^/dC^  vs  ) .  although  they 
represent  configurations  markedly  different  to  the  civil  transport,  the  results  from 
which  are  shown  in  Figs  2  to  5. 

5  CONCLUSIONS 

A  rapid  and  reliable  method  for  estimating  equivalent  parabolic  drag  polars  from 
wind-tunnel  data  has  been  described,  and  compared  with  an  alternative  technique  described 
in  Ref  3.  The  new  method  has  some  advantages  over  that  in  Ref  3,  in  that  the  part  of  the 
measured  drag  polar  that  may  be  represented  accurately  by  a  parabola  is  identified  at  the 
outset.  An  important  part  of  the  current  method  is  the  use  of  an  interactive  computer 
graphics  system,  where  linear  portions  of  the  polar  slopes  may  be  selected  from  the  key¬ 
board  for  processing  by  the  computer. 

If  such  a  graphics  system  is  available  in  the  wind-tunnel  control  room,  then  the 
possibility  arises  of  performing  the  separated-f low  wake  blockage  correction  'online', 
i.e  as  data  is  acquired.  For  example,  in  a  typical  test  to  determine  the  maximum  lift 
coefficient  of  a  model,  the  test  proceeds  by  measuring  forces  over  a  range  of  incidences, 
increasing  incidence  between  successive  datapoints.  If  dC^/dC^  is  plotted  against 
as  the  test  proceeds,  the  test  controller  can  observe  when  a  linear  portion  of  the  curve 
appears,  and  signal  the  computer  to  use  this  linear  portion  to  fit  a  parabolic  drag  polar. 
The  test  results  can  then  be  corrected  for  wake  blockage  effects  due  to  separated  flow. 
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■'  A  method  for  obtaining  a  'best-fit'  parabolic  drag  polar  from  measured  force 
results  is  presented.  Tbe  method  is  used  in  estimating  the  blockage  correction  due 
to  separated  flow  that  must  be  applied  to  wind-tunnel  results. 

An  important  part  of  the  method  is  the  use  of  computer  graphics  to  identify 
the  part  of  the  measured  drag  polar  that  may  best  be  described  by  a  parabola, 
and  thus  used  as  the  basis  for  a  curve  fit. 

The  method  is  fast,  when  co^ared  to  existing  curve-fitting  techniques, 
and  suited  to  use  with  an  ^nline^data  reduction  system. 


